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4^ Abstract: 



X 



The internal spin structure of the A is of special importance for the under- 



$_i ■ standing of the spin structure of hadrons in general. The comparison between 



the nucleon and A allows for a test of the relevant flavour-symmetry break- 
ing effects. Using nonperturbatively 0(a) improved Wilson fermions in the 
quenched approximation we have calculated the first moments of the unpo- 
larised, longitudinally polarised and transversity quark distribution functions 
in the A. The results indicate that flavour symmetry breaking has little effect 
on the internal spin structure, in accordance with model based expectations. 
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1 INTRODUCTION 



The detailed investigation of the nucleon spin structure during the last few 
years has fuelled a large number of efforts to understand the results with 
quark models of various kinds. It turned out during the course of these in- 
vestigations that SU(3) flavour breaking is of central importance for any such 
effort. A small selection of relevant articles can be found in [|l|. Presently 
the different predictions cover a rather broad range, and while some authors 
claim great accuracy, others conclude that the existing experimental data on 
hyperon decays is insufficient for any stringent prediction ||. The A spin 
structure, in turn, is especially sensitive to flavour SU(3) breaking R |1 H. 
While in the naive SU(6) quark model the spin of the A is carried exclusively 
by the s quarks, the SU(3) rotated results for the nucleon spin structure 
suggest that the s and s quarks carry only rs 60% of the A spin while u, u, 
d and d quarks quarks contribute ~ —40%. Any specific assumptions about 
SU(3) symmetry breaking will strongly affect these conclusions. Obviously 
additional input is needed, some of which we want to provide with the lattice 
results presented in this paper. 

The spin structure of the A is also experimentally of special interest, be- 
cause its polarisation can easily be measured via the self-analysing weak 
decay A — > pn~ . Indeed, the A polarisation has been determined at the 
Z° pole in e + e~ annihilation || where it is mainly due to the s quarks, 
which when produced via Z° decays have an average polarisation of —0.91. 
Furthermore, A polarisation has been studied in deep-inelastic scattering of 
polarised positrons on unpolarised protons in the current fragmentation re- 
gion, i.e. selecting As which most likely originate from the struck quark 0. 
Needless to say, however, that the interpretation of the experiments is com- 
plicated by fragmentation effects, As from the decay of heavier hyperons, etc. 
A and A polarisation was also analysed for charged- current neutrino nucleon 
reactions ||. A completely unsettled puzzle is the very strong polarisation 
of As produced in unpolarised p + N reactions || . The knowledge of the 
internal spin structure of the A is certainly a necessary ingredient for any 
explanation. 

Thus, there is ample motivation to perform lattice calculations which provide 
information on the internal A structure. For a spin 1/2 baryon the fraction 
Aq of the spin carried by the quarks (and antiquarks) of flavour q is given in 
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terms of the forward matrix element of the axial vector current: 

(P, s|9%75#, s ) = 2s M Ag , (1) 

where s 2 = —M 2 (M = mass of the baryon) and the states are normalised 
according to (p,s\p',s') = 2E p (2n) 3 5(p — p')5 ss >. Simultaneously we have 
calculated the transversity matrix element (the tensor charges Sq): 

2 

(p, s\q ia/tujs q\p, s) = -T7{ s ^Vu ~ s u p^)5q . (2) 

The physical origin of the difference between Aq and 5q is that in the first 
case the baryon is boosted along its spin direction and in the second case 
perpendicular to it. Thus this difference has a very specific sensitivity to 
details of the internal baryon structure, being e.g. zero in the non-relativistic 
limit. Operators with additional derivatives can be used to extract also some 
higher moments of the corresponding distribution functions. 
Our aim is to study how strongly the mass difference between the strange 
quark and the light quarks affects the flavour SU(3) symmetry between the 
(polarised) distribution functions in the nucleons and the A. Perfect SU(3) 
symmetry predicts the equations 

As A = (2Au p - Ad p + 2 As p ) /3 , 
Ad A = Au A = (Au p + 4Ad p + As p )/6 (3) 

and the same relations for the tensor charges. 

A preliminary account of some of our results has already been given in Ref. 

E3- 



2 THE SIMULATION 

We have performed quenched simulations with the Wilson gauge action at 
(3 = 6.0 using nonperturbatively 0(a) improved fermions (clover fermions) 
with csw = 1-769. The lattice size was 16 3 x 32. We worked with nine 
combinations of hopping parameters: k u = k s e {0.1324,0.1333,0.1342} 
corresponding to bare quark masses of ~ 166, 112, 58 MeV, respectively. So 
we can extrapolate to the chiral limit in k u = and interpolate in k s to the 
physical value k*. For the critical hopping parameter k c we take the value 



0.135201 determined from the PCAC quark mass [11|, k* = 0.1341 was fixed 
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by requiring the pseudoscalar mass m PS for k u = n d = k c and k s = k* to be 
equal to the K + mass of 494 MeV (with the scale set by the force parameter 
tq = 0.5 fm). In the extra- and interpolation we assumed a linear dependence 
of mp S on the quark mass. 

As an interpolating field for the A we used (employing Euclidean notation 
from now on) 

&a(t) = E e ijkSia(x) (uj (x)C^f 5 d k (xfj (4) 

X, X4=t 

where C is the charge conjugation matrix; i, j, and k are colour indices and 
a is a Dirac index. The quark fields were (Jacobi) smeared to improve the 
overlap with the A state. 

As usual, the bare matrix elements are determined from ratios of three-point 
functions and two-point functions 

with appropriate spin projection matrices T and 17 The choice for V and T 
as well as t and the smearing parameters can be found in ||12|| . The quality 



of the data can be judged by how well pronounced a plateau is obtained for 
this ratio as function of r. 

On the lattice the choice of the operators O is a non-trivial task, because the 
discretisation reduces the symmetry group of (Euclidean) space-time from 
0(4) to the hypercubic group H(4) C 0(4). Therefore, one has to find 
combinations of operators which avoid the problem of operator mixing JT3| . 



After renormalisation, the matrix elements of these operators are expressed 
in terms of the reduced matrix elements t>2, cto, a i and to which correspond to 
moments of the parton distribution functions. We work with the operators 

Ov 2 = \q (74 Da ~\ (li Si +72 B 2 +73 D3X) q , (6) 
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O ao = qi2l5q , (7) 

O ai = ^ q (^7475 D 2 +7275 D*j q , (8) 

O to = qa 2 4l5q ■ (9) 
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While the w 2 operator Eq.(|GJ) (the v 2 b operator in the notation of [0) is used 
in a spin-averaged matrix element, the other operators require a polarised 
baryon state to give a non-zero result. In these latter cases we choose the 
spin to point into the 2-direction. In parton model language we have: 

A = (xq(x)) , a q Q = 2(Aq(x)) , a\ = 2(xAq(x)) , ig = 2(5q(x)) . (10) 

In order to reduce the cut-off effects from 0(a) to 0(a 2 ) also in matrix 
elements, the improvement of the fermionic action has to be accompanied by 
the improvement of the operator under study. For the axial vector current 
the improved renormalised operator has the form 

A <mpr = Zao(1 + baoam )(A, + ac ao d,P) (11) 

with the bare axial vector current A^(x) = q(x)jfj, r Y5q(x) , the pseudoscalar 
density P(x) = q^^q^x), and the bare quark mass am = 1/(2k) — 1/ (2k c ). 
The improvement term d^P vanishes in forward matrix elements, so we do 
not need the coefficient c ao . For the other operators one can use the equation 
of motion to define the improvement coefficients analogous to c aQ such that 
they are of order g 2 and and thus small. We therefore decided in the end 
to neglect the improvement terms altogether. For baryons with momentum 
zero the ratios Rq are related to the reduced matrix elements (|K]) by 

R V2 — 7 ; r — Mt>2 , 

Z V2 (\ + b V2 am) 2k 



Rn 



i 1 1 



Z ao (l + b ao am) 2k 2 

i 1 M 

Rai ~ ~Z ai (l + b ai am)2K~T ai ' 
111 
Z to (l + b to am) 2k 2 

Depending on the flavour of the quark one has to insert for k either k^ or 
k s . These relations account for the operator renormalisation and different 



normalisations on the lattice and in the continuum |12| [14| . The renormal- 
isation constants Z V2 , Z ao , Z ai , Z to depend on the renormalisation scale /i 
of the continuum theory and the lattice cut-off 1/a or equivalently f3. One 
exception is Z ao , which depends only on /3, because in the continuum the 
anomalous dimension is zero due to current conservation. The scale depen- 
dence can be factorized into a factor depending on /i and one depending on f3. 
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For the renormalisation constants we employ a variation of tadpole improved 
perturbation theory - TRB-PT - as briefly described in |15|]. The renormali- 
sation group invariant form is first found (accurate to two loop perturbation 
theory) and then converted back to the MS-scheme at /i 2 =4 GeV 2 , using 
two-loop perturbation theory for t an d three-loop perturbation theory for 



ai and f 2 . We obtain |T6] 



Z V2 (fx 2 


= 4 GeV 2 




6.0) = 


1.11106 


K 


= 1.25803 






w= 


6.0) = 


0.83232 


K 


= 1.27134 




= 4 GeV 2 


,p = 


6.0) = 


1.11800 


b ai 


= 1.24313 


z uAv 2 


= 4 GeV 2 


,p = 


6.0) = 


0.88924 


K 


= 1.24626 



(13) 

Z ao and b ao have been determined non-perturbatively (on the lattice) in (T7J 
with the result 

Z ao (P = 6.0) = 0.807(2) (8) b ao = 1.28(3) (4) . (14) 



This value for Z ao deviates by « 4% from the perturbative value in Eq. ( |I~3"D . 
We used for a Eq.(|T4"D and otherwise Eq.(^). At this point we also want 
to address the question with which phenomenological distribution functions 
we should compare. Our Z values include all contributions of order a s and 
some (but not all) contributions of higher order. This suggests that we should 
compare them with an NLO (or higher) DGLAP-fit to the experimental data 
in the MS scheme and for the scale 4 GeV 2 . 



3 RESULTS 

From the A masses at our nine combinations of k,j, k, s we have computed A 
masses at = k c by linear extrapolation of M\ in 1/tZd- These 12 masses 
are plotted in Fig. [I]. The filled triangles denote the chiral limit k^, k u — > 
k c for each k s . Our value for k*, which we fixed by M K , reproduces quite 
accurately the ratio M\/M p = 1.19. This fits in nicely with the observation 
that Ms/M p and M^/M p are also rather well reproduced by quenched sim- 
ulations (see, e.g., JTTJ] ) . The plot shows clearly the breaking of the SU(3) 
flavour symmetry. In particular, the dependence of the masses on Kd is rather 
pronounced. 

In Fig. we plot our bare results, i.e. the results without the factors 



Z(l + bam) in Eq. (12), for (Aq(x)), (xAq(x)) (xq(x)), and (Sq(x)) in the 
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04 3.0 




7.3 7.35 7.4 7.45 



Figure 1: The square of the ratio M\/M p with the proton mass M p taken in 
the chiral limit versus 1/k s . The different symbols correspond to the different 
values of including the chiral limit. The crosses (left to right) indicate the 
physical value of M A /M p , 1/k c , and respectively. 
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u, d 


s 


(xq(x)) 


0.20(1) 


0.27(1) 


(Aq(x)) 


-0.02(4) 


0.68(4) 


(xAq(x)) 


0.033(8) 


0.23(2) 


(5q(x)) 


0.03(2) 


0.74(4) 



Table 1: Our results for different moments of polarised and unpolarised 
quark distribution functions in the MS scheme at /i 2 = 4 GeV 2 . 



A versus l/« s (q = s, d; the results for u are the same as for d). We plotted 
the bare results because these are the numbers we actually used for inter- 
and extrapolation. In contrast with the case of M\/M p , the dependence on 
Kd is rather weak. The values corresponding to Kd = k c have been obtained 
by extrapolating the bare matrix elements linearly in (filled triangles). 
Finally, by interpolating the bare matrix elements for Kd = k c linearly in 
1/k s to 1/fc* we obtain the desired A matrix elements. The numbers plotted 
in the figure are given as tables in the appendix. The physical results for 
the renormalisation scale fi 2 = 4 GeV 2 are collected in table |l|. Thus we find 
that the d and u quarks carry a somewhat smaller momentum fraction than 
the s quarks, 20% each as compared to 27% which is quite intuitive in view 
of the larger s mass. We also find that the tensor charge Sq is close to Aq for 
the s quark, which again fits to our understanding of transversity, as both 
should be identical in the non-relativistic limit. Note that a positive value 
of (xAd(x)) and a smaller or even negative value for (Ad(x)) imply a sign 
change of Ad(x) as function of x. This would agree with model predictions 

These numbers are subject to systematic uncertainties which we will now 
discuss. For the nucleon similar quenched lattice calculations always gave 
quark momentum fractions which were too large by typically about 20% 
and this disagreement was usually attributed to quenching and the neglect 
of quark-line disconnected contributions. Recently it has been pointed out 
that the extrapolation to the chiral limit could be of similar or even greater 
importance [frS|| . If this is true we expect that the corrections from a nontriv- 



ial chiral extrapolation can affect substantially the values for Aq. Another 
source of uncertainty is the lacking continuum extrapolation as the current 
calculations have been performed at one lattice spacing only. 
Let us now discuss our results for Aq in more detail on the basis of table ||. 
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Figure 2: The bare results for (Aq(x)), (xAq(x)), (xq(x)), and (5q(x)) in 
the A versus 1/k s . The results for the different values of n d ^ k c have been 
slightly displaced horizontally. 
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Aua = Ag?a 


As A 


quark model 





1 


exp. + SU(3) F 


-0.17(3) 


0.63(3) 


MC + SU(3) F 


-0.016(9) 


0.65(2) 


this work 


-0.02(4) 


0.68(4) 



Table 2: Comparison of our results for the longitudinal quark polarisation 
in the A with the naive quark model prediction, the flavour SU(3) rotated 
experimental values for the proton and flavour rotated lattice results for the 
proton. 

In this table we give in the third line the A matrix elements as they follow 
from our Monte Carlo results for the proton matrix elements by the use of 
SU(3)f (see Eq.(3)). They agree quite well with the matrix elements com- 
puted directly (fourth line). This implies that the flavour symmetry breaking 
effects in the matrix elements are rather small, which can also be concluded 
from Fig. |2|. As the mass difference between the light quarks and the strange 
quark was consistently taken into account this provides a strong argument 
that the A and proton spin structures are, in good approximation, simply 
related by an SU(3) transformation. Consequently, the values given in the 
second line, which were computed from the proton spin structure under the 
assumption of flavour SU(3)(see, e.g., ||), should be quite reliable. 
Because we performed a quenched calculation it might be more consistent 
to compare our results to the prediction for the valence quark contribution, 
e.g. by Ashery and Lipkin ||. They obtained Au\ = Ad\ = —0.07(4), 
Asa = 0.73(4). Notice that all results differ markedly from the predictions 
of the (naive) quark model shown in the first line. 

4 SUMMARY 

We have studied various aspects of the A (spin) structure taking the mass 
difference between the quarks into account. We found that our results agree 
nicely with general expectations: The momentum fraction carried by the s 
quark is larger than for light quarks, the values for the tensor charge 5s and 
As are similar. Our main result is, however, that SU(3) flavour symmetry 
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appears to be much less violated in the matrix elements than in the baryon 
masses, which is in agreement with the empirical observation that the hy- 
peron semileptonic decays can be parametrised rather well assuming SU(3) 
flavour symmetry. Presently the use of As as a probe of the nucleon spin 
structure in e.g. deep-inelastic lepton-nucleon scattering still has to rely on 
models, but the properties of the latter should definitely be checked by a 
comparison with our lattice results. 

We expect that our results have a similar uncertainty as those for the (spin) 
structure of the proton and that the origin of this uncertainty is the same. 
When configurations with lighter dynamical quarks become available and our 
understanding of chiral extrapolation improves, high precision results for the 
A should become reachable along the lines discussed in this contribution. 
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APPENDIX 

The following tables contain the numerical values for our lattice results. 
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af a\ 


0.1324 0.1324 
0.1333 0.1324 
0.1342 0.1324 
0.1324 0.1333 
0.1333 0.1333 
0.1342 0.1333 
0.1324 0.1342 
0.1333 0.1342 
0.1342 0.1342 


-0.02(1) 1.64(2) 
-0.02(2) 1.65(3) 
-0.05(4) 1.66(4) 
-0.01(2) 1.61(3) 
-0.02(3) 1.63(3) 
-0.05(5) 1.63(5) 
-0.00(3) 1.58(4) 
-0.01(5) 1.59(6) 
-0.05(7) 1.60(9) 


0.014(4) 0.350(7) 
0.019(5) 0.378(9) 
0.027(7) 0.416(15) 
0.017(4) 0.321(8) 
0.024(6) 0.350(10) 
0.037(9) 0.387(18) 
0.022(6) 0.291(9) 
0.032(7) 0.318(13) 
0.051(12) 0.357(27) 


k c 0.1324 
k c 0.1333 
k c 0.1342 


-0.05(4) 1.67(5) 
-0.05(6) 1.64(6) 
-0.06(10) 1.61(10) 


0.03(1) 0.45(2) 
0.05(1) 0.42(2) 
0.06(2) 0.38(3) 


K c K g 


-0.06(9) 1.62(9) 


0.06(2) 0.39(3) 



Table 3: Measured bare values for a and a\. 





Kg 


v d 2 


V 2 


l o 


l o 


0.1324 


0.1324 


0.188(2) 


0.207(3) 


0.02(1) 


1.58(2) 


0.1333 


0.1324 


0.179(3) 


0.223(4) 


0.02(1) 


1.60(3) 


0.1342 


0.1324 


0.169(5) 


0.244(5) 


0.04(2) 


1.61(3) 


0.1324 


0.1333 


0.194(3) 


0.193(4) 


0.03(1) 


1.57(3) 


0.1333 


0.1333 


0.186(4) 


0.209(4) 


0.03(1) 


1.58(3) 


0.1342 


0.1333 


0.177(7) 


0.230(6) 


0.05(3) 


1.59(4) 


0.1324 


0.1342 


0.202(3) 


0.177(5) 


0.03(1) 


1.55(5) 


0.1333 


0.1342 


0.195(4) 


0.192(6) 


0.04(2) 


1.57(6) 


0.1342 


0.1342 


0.188(9) 


0.213(9) 


0.08(4) 


1.57(8) 


K c 


0.1324 


0.159(7) 


0.262(6) 


0.04(3) 


1.62(4) 


K c 


0.1333 


0.168(8) 


0.248(7) 


0.06(4) 


1.61(5) 


K c 


0.1342 


0.180(11) 


0.229(10) 


0.08(5) 


1.59(10) 


K c 




0.178(10) 


0.233(10) 


0.08(5) 


1.59(8) 



Table 4: Measured bare values for v 2 and to. 
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